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ABSTRACT 

on" 

Qv^ I In the external shock model, gamma-ray burst (GRB) emissions are produced by 

the energization and deceleration of a thin relativistic blast wave due to its interactions 
1-^ ■ with the circumburst medium (CBM). We study the physical properties of an analytic 

Ft ■ function which describes temporally-evolving GRB spectra in the limit of a smooth 

CBM with density n{x) oc x~'^, where x is the radial coordinate. The hard-to-soft 
Cn , spectral evolution and the intensity-hardness correlation of GRB peaks are reproduced. 

We predict that (1) GRB peaks are aligned at high photon energies and lag at low 
K* I energies according to a simple rule; that (2) temporal indices at the leading edge of a 

Q I GRB peak display a well-defined shift with photon energy; and that (3) the change 

rn ' in the spectral index values between the leading and trailing edges of a GRB peak 

f— ^ ■ decreases at higher photon energies. The reason that GRBs are usually detected with 

^N . vFi, peaks in the 50 keV - several MeV range for detectors which trigger on peak flux 

over a fixed time interval is shown to be a consequence of the inverse correlation of 
peak flux and duration of the radiation emitted by decelerating blast waves. 
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c/3 . Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanisms: 

. . , nonthermal 
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Cd ' 1. Introduction 

The Beppo-SAX results have revolutionized our understanding of GRBs by opening a window 
on X-ray, optical, and radio afterglows which are delayed from the prompt X-ray and 7-ray 
emissions by several hours and more (e.g., Costa et al. 1997; Feroci et al. 1998; Piro et al. 1998; 
van Paradijs et al. 1997; Frail 1998). The decaying power-law long- wavelength afterglows are 
compellingly explained as the emissions from a relativistic blast wave which decelerates and 
radiates through the process of sweeping up material from a uniform CBM (e.g., Paczyiiski &; 
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Rhoads 1993; Katz 1994; Meszaros & Rees 1997; Wijers, Rees, k Meszaros 1997; Tavani 1997). 
Afterglow behaviors more complicated than simple power laws may involve effects of extinction, 
self-absorption, scintillation, and CBM structure. 

Meszaros & Rees (1993) proposed that the energy of a fireball's blast wave could be efficiently 
converted into radiation during the prompt gamma-ray luminous phase of a GRB through the 
interaction of a blast wave with the CBM. In the external shock model, all GRB emissions are 
due to the effects of a single thin blast wave which interacts with the CBM. Because the directed 
kinetic energy of the blast wave is converted into nonthermal particle energy by sweeping up 
material from the CBM, the emitted radiation depends crucially on the density distribution 
of the CBM in the blast wave's path. Smoother GRB time profiles therefore represent more 
uniform CBMs, at least within the Doppler cone from which most of the detected GRB emissions 
originate. Conversely, erratic and spiky GRB time profiles represent blast-wave deceleration in 
highly textured CBMs (Dermer & Mitman 1999). 

In a recent paper (Dermer, Chiang, & Bottcher 1999; hereafter DCB), we proposed a 
parametric description of the radiation observed from a spherically expanding blast wave which 
decelerates and is energized by sweeping up material from a smooth CBM. By "smooth" we mean 
that the CBM density distribution is adequately represented by the expression n{x) = riQX^^, 
where x is the radial coordinate. In this limit, the deceleration of the blast wave produces a time 
profile which mimics the so-called Fast Rise, Exponential Decay (FRED) GRB light curves. Our 
ansatz is that a smooth CBM produces an ideal FRED light curve. Conversely, a GRB which 
exhibits the classic FRED-type profile results from a fireball embedded within and expanding into 
a smooth CBM. 

In this Letter, we make three quantitative predictions that can be tested using BATSE and 
Beppo-SAX data from bright FRED-type GRB time profiles, though the predictions would be 
most thoroughly tested with a GRB telescope which is much more sensitive than BATSE and 
measures prompt GRB emission in the range between ~ 1 keV and several MeV. In Section 2, 
we describe the model equations derived in DCB. Spectral and temporal predictions are given 
analytically and described graphically for a model fireball with Fq = 300. In Section 3 we show 
that the reason that GRBs are usually observed with i^Fu peaks in the ^ 50 keV-several MeV 
range is understood when account is taken of blast-wave physics and the triggering properties of 
burst detectors. 



2. Analytic Description and Predictions of the Blast Wave Model 

In DBC, a parametrization of the nonthermal synchrotron emission from a decelerating blast 
wave observed at time t and photon energy e = hu/nieC^ was proposed. The function is based on 
the analytic approach of Dermer & Chiang (1998) and the numerical results of Chiang & Dermer 
(1999), which were in turn based upon the blast-wave physics developed by, e.g., Blandford & 



McKee (1976), Rees & Meszaros (1992), Piran & Shemi (1993), Meszaros & Rees (1993), Meszaros, 
Laguna, & Rees (1993), Waxman (1997), and Vietri (1997). 

The deceleration of the relativistic blast wave is assumed to follow the expression T{x) = Fq 
when X < Xd, and r(a;) = TQ{x/xd)~^ when x^ < x < x^^q . Here Fq is the initial bulk Lorentz 
factor of the blast wave and the decleration radius Xd = 2.6 x 10"*^^ [(1 — r7/3)£^54/(n2F3QQ)]^'^ 
cm. In this relation, the burst source emits dE/dQ = W^'^E^^/^Att) ergs sr~-^, Fq = 3OOF300, and 
1^2 = ^0/(10^ cm~^). For simplicity, we consider only spherically symmetric blast waves. The 
parameter g specifies the radiative regime, and g ^ {3 — r])/2 and g ^ 3 — ij in the adiabatic and 
radiative limits, respectively. The deceleration time scale t^ = 9.7(1 + z)[{l — 7?/3)£'54/(n2F3QQ)]^/^ 
s. 

The function proposed in DBC to model the blast-wave radiation has the broken-power law 
form 

where v and 6 are the lyF^ spectral indices at energies below and above the temporally-evolving 
break energy ep(t), respectively. Expressions for ep{t) = £()[T{x) /Tq]'^{x /xd)~^'^ and the 
time-varying amplitude Pp{t) = IlQ[T{x)/TQ]^{x/xd)'^^^ are given in DBC by eqs.(9) and (16), 
respectively. Eq. (1) depends on the specification of the nine parameters listed in Table 1. The 
equipartion parameter q = [Ch{'''/^)]^ Ce^ where r is the compression ratio and ^h and ^e are the 
magnetic field and electron equipartition values, respectively (see eqs. [7] and [8] in DCB). The 
parameters v, 6, g, ij, and q are assumed to be independent of time. This assumption is obviously 
not true in general, but is realized in the limit that high-quality broadband measurements 
from a bright FRED-type GRB are made over a sufficiently short observing times 5t so that 
radiative-cooling and magnetic field-evolution time scales in the comoving frame are ^ T5t/{l + z). 

The thick curves in Fig. 1 show time profiles calculated at various observing energies e, using 
the equations of the previous section. The thin curves show spectral indices calculated between 
e and 2e. Here we show results for fireballs with Fq = 300 which are located at redshift z = 1; 
other parameters are listed in Table 1. In this example, td = 9.6(1 + z) s and Sq = 2.43/(1 + z). 
The model GRB light curves shown in Fig. 1 display a rapid rise followed by a gentler decay 
which approaches a power-law afterglow behavior oc t^^-^"^ (see eq. [22] in DCB) at late times 
t ^ td- The overall shape of the model light curves shown in Fig. 1 resembles the characteristic 
FRED-type GRB light curve (e.g., Fishman & Meegan 1995). The model GRB represents a very 
bright BATSE GRB if located at z = 1, but would fall below BATSE detectablity if located at 

The analytic approximation used for the evolution of F(x) gives spectral index curves which 
are constant when t < td- At later times, the blast-wave deceleration means that an observer 
measuring a photon spectrum over a fixed range of energies will sample photons which are 
produced at progressively higher energies in the blast wave's frame of reference. Because a 
nonthermal synchrotron spectrum from a power-law distribution of electrons with a low energy 



cutoff is very hard at low energies and softens at higher energies, blast-wave deceleration produces 
the hard-to-soft spectral evolution observed in many GRB time profiles (Norris et al. 1986). If the 
observer is monitoring a GRB at e ;^ £q, the flux is brightest when t ~ t^. As can be seen from 
Fig. 1, the spectrum remains hard until t ';^ td, after which the blast wave begins to decelerate 
and the spectrum softens. Thus the spectrum is hardest when it is most intense, accounting for 
the hardness- intensity correlation observed in GRB light curves (Golenetskii et al. 1983). If the 
blast wave produces short time scale variability by interacting with inhomogeneities in the CBM 
(Dermer &: Mitman 1999), then the individual pulses in GRB profiles would likewise exhibit 
spectral hardening and subsequent softening, as generally observed in well-defined pulses of GRB 
light curves (Ford et al. 1995; Crider et al. 1998). The qualitative ability of the blast wave model 
to explain these empirical trends has been pointed out previously by Panaitescu & Meszaros 
(1998). 

We now propose three quantitative predictions to test the validity of this model. As is evident 
from Fig. 1, the peak flux shifts to later times at lower photon energies. The peak times are given 
analytically by the expression 
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(see also Chiang 1998). By plotting the times of a well-defined peak for a GRB measured over a 
large energy range, one tests the external shock model and constrains values of Sq and the index 
{2g + l)/{4g + ri/2). The peak shifting will be more pronounced in GRBs with larger values of £q. 

Our second prediction is that the temporal indices of a pulse profile vary with photon energy 
according to the relation xi^'ii) = dlnP{e,t)/dlnt. This yields a complicated analytic expression 
from eq. (1), though it is easily calculated numerically. When t < td, however, the analytic form 
gives 

f (1 + lMr^it/tdf'^n+vv/2 for e « £0, and 

P(e,t) = no X I (3) 

[ (1 + |)(i)"'(Vid)'"''""'/' for e » £0. 

Thus we predict that the temporal index x changes from 2 — 7/(1 — v/2) at low energies to 
2 — 77(1 -|- (5/2) at high energies. For a uniform CBM, there is no shift and the temporal index 
X = 2. When rj j^ 0, the relation between the temporal and spectral indices indicated by eq. (3) 
implies the value of rj. 

As can also be seen from Fig. 1, the spectral index is hardest at the leading edge of a GRB 
peak and softens at the trailing edge, with the change in the values of the spectral index decreasing 
towards higher photon energies. This prediction can be made quantitative by taking the derivative 
of eq. (1), giving the result 
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for the vFi, spectral index, where y = e/ep{t). The external shock model can be tested by plotting 
the energy-dependent variation of the spectral index across the peak of a GRB from eq. (4). 

Although the above equations provide a simple analytic characterization of the predictions of 
the external shock model, they depend crucially on the r(x)-prescription noted above. Moreover, 
electron cooling is not taken into account in the analytic model. We used a numerical simulation 
code (Chiang & Dermer 1999) to determine the reliability of the analytic predictions. The 
results are shown in Fig. 2 using the parameters of Table 1. Because the radiative regime g in 
the simulation depends on the magnetic field H{G), the value of £^h = -ff^/[327r?72pC^no(r/4)] 
was adjusted until the bulk Lorentz factor approached the asymptotic behavior T oc x~^ in the 
deceleration phase. This occurred for ^h = ?, x 10^^. The injection spectrum of electrons is 
chosen oc ^^•^•^ in order to give an uncooled spectrum with 6 = 0.2. Electron cooling is seen to be 
important at high photon energies. The peak of the light curves at large photon energies occurs 
at ~ id/2, with a peak flux ~ 2x greater than the analytic estimate. We find that the qualitative 
trends of the predictions described above are reproduced in the detailed calculation. The analytic 
model can be used to confine the parameter range, but detailed fitting to data should employ the 
more accurate numerical simulations. 



3. Resolution of the Relativistic Beaming/i/Fi^ Peak Energy Paradox 

Even if cosmic fireballs were produced with a narrow range of values of Fq, their uFi^ peaks 
would be distributed over a wide range of photon energies since Sq oc Eg. It has therefore been 
something of a mystery to understand why GRBs are usually detected with uFi, peaks in a narrow 
range of photon energies between ~ 50 keV and ~ 1 MeV (see, e.g., Malozzi et al. 1995; Piran 
& Narayan 1996) without requiring fine-tuning of the parameters. Brainerd's (1994) Compton 
attenuation model, for example, has been advanced to specifically address this puzzle. By contrast 
to GRBs, the nonthermal emissions in blazars are also thought to originate from relativistic 
outfiows, yet they have vFy peaks in both their synchrotron and Compton components which 
range over three orders-of-magnitude or more in photon energy (e.g., von Montigny et al. 1995). 

Using the spectral form proposed in DCB, we show how this paradox is resolved for a uniform 
CBM (the generalization to 77 7^ is straightforward). Suppose an instrument triggers on peak 
flux measured over the time scale At and over a narrow range of photon energies centered at 
€(1- When the mean duration tp{ed) of the fireball measured at e^ is longer than At, then the 
detector triggers on peak fiux. Erom the expressions for ep{t) and Pp{t), we find that fireballs 
with Eq = Eq — 240 [(1 + z)€dl{n2 9-3)]^''^ are observed with the peak of their vF^, spectrum 
at photon energy e^ at the moment when the received bolometric power is greatest (i.e., when 
t = td). Dirtier fireballs with Eq < Eq produce a peak flux at ed which declines according to 
the relation Pp{ed,tp) oc Eq , as can be derived by inserting eq. (2) into eq. (1) in the limit 

^d ^ £0 (see also eq. [19] in DCB). Because 6 > 0, the peak fiux therefore rapidly decreases with 



decreasing Fq. Cleaner fireballs with Fq > Fq are also more difficult to detect because in this case, 
Ppi^d^tp) oc Fq^ ' (cf. eq. [20] in DCB). In order to be observed at all, blast waves cannot be 
perfectly adiabatic, so that g > 3/2 and the peak flux detected with a GRB instrument decreases 
with increasing Fq. 

Cleaner fireballs are more difficult to detect than indicated above, because when td <C At, 
the detector triggers on fluence rather than flux. The fluence F decreases with increasing Fq 
according F oc tpPp{eii,tp) oc Fq x Fq^ , using eq. (2) and the expression for ep{t) with 

r/ = 0. Because of the emission properties of blast waves, detectors which trigger on peak flux 
over a fixed time window will therefore be most sensitive to fireballs which have vFy peaks in the 
energy range of the detector. Hence the discovery of new classes of clean and dirty fireballs must 
consider varying time windows and energy ranges in the triggering criteria of a detector, and must 
also contend with the levels of background radiation as discussed in DCB. 

In separate work (Bottcher &: Dermer 1999), we model the triggering properties of GRB 
detectors and show that compared to the BATSE GRB trigger rate, the dirty fireball rate is 
poorly known due to selection biases against their detection. Observational analyses of Ariel V 
X-ray data (Grindlay 1999; see discussion in DCB) provides the strongest available limit on the 
frequency of dirty fireballs. By contrast, the relative rate of clean fireballs is strongly constrained 
in the analysis of Bottcher &: Dermer (1999) unless there is a strong anti-correlation between Fq 
and q. 

In summary, we have used a simple analytic description of the external shock model of GRBs 
to pose three predictions which can be tested by fitting eqs. (2)-(4) to high quality data obtained 
with a GRB telescope sensitive in the ~ 1 keV - several MeV range. Although we have employed a 
simple mathematical parametrization, detailed fits to data should use numerical simulations of the 
decelerating blast wave (e.g., Chiang & Dermer 1999, Panaitescu &: Meszaros 1998). We have also 
explained why detectors which trigger on peak flux are most sensitive to fireballs which produce 
I'Fi, peaks in the 50 keV - several MeV range, resolving a long-standing puzzle in relativistic 
beaming models of GRBs. 
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Fig. 1. — Light curves (thick hnes) and spectral indices (thin hnes) from a model GRB produced 
by a fireball with Fq = 300 located at redshift z = 1. Other parameters are listed in Table 1. 
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Fig. 2. — Numerical simulation results corresponding to the analytic model shown in Fig. 1 (see 
text for details). Equipartition between the energy in electrons, protons, and magnetic field is 
assumed. 



